Complex Periodic Orbits and Tunnelling in Chaotic Potentials 
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We derive a trace formula for the splitting-weighted density of states suitable for chaotic potentials 
with isolated symmetric wells. This formula is based on complex orbits which tunnel through 
classically forbidden barriers. The theory is applicable whenever the tunnelling is dominated by 
isolated orbits, a situation which applies to chaotic systems but also to certain near-integrable ones. 
It is used to analyse a specific two-dimensional potential with chaotic dynamics. Mean behaviour 
of the splittings is predicted by an orbit with imaginary action. Oscillations around this mean are 
obtained from a collection of related orbits whose actions have nonzero real part. 

PACS numbers: 03.65.Sq, 73.40Gk, 05.45. +b, 0.320. +i 



In this letter we develop a method for computing tun- 
nelling effects in quantum states associated with chaotic 
regions of phase space. While the influence of chaotic 
dynamics has been actively studied in the context of 
chaos-assisted tunnelling between EBK-quantised tori 
|0,D, much less attention has been paid to the problem 
of tunnelling between chaotic states. This is presumably 
because there is at present no semiclassical theory for 
individual chaotic eigenstates. However, there does ex- 
ist a well-developed theory for the spectral properties of 
such systems in terms of periodic orbits jsj . We will show 
how this can be extended to a calculation of spectral tun- 
nelling averages using complex periodic orbits. 

Complex multi-dimensional trajectories were intro- 
duced in Q to understand barrier penetration problems 
as motivated by quantum field theories. Complex peri- 
odic orbits were first used by Miller in one-dimensional 
potentials, to derive splittings and resonance widths and 
these calculations were extended in |^ to the rotational 
spectrum of SFg. In the context of chaotic maps, com- 
plex periodic orbits were used to calculate band gaps . 
Complex trajectories in chaotic maps have also been ex- 
plored in the time domain Q . 

A commonly used probe for tunnelling effects is to 
study the spectra of two symmetric wells which have a 
barrier between them. When classical trajectories are 
localised within one well or the other, one finds that en- 
ergy levels come in symmetric-antisymmetric pairs E^, 
with small splittings between them. Denote the split- 
tings by AEn and the mean levels by En, so that E^ — 
En T AEn/2. A standard periodic orbit calculation us- 
ing real orbits yields a set of doubly-degenerate levels 
En , approximating the mean levels £"„ . The degeneracy 
arises because each periodic orbit has a symmetric part- 
ner in the opposing well. Therefore, the most naive use 
of periodic orbit theory fails to predict the existence of 
splittings. However, since including complex orbits does 
yield splittings in one dimension |^,^, we are motivated 
to pursue this approach in higher dimensions. 

A literal extension of this work to chaotic problems is 
too difficult because it involves finding tiny differences 



between poles in traces or zeros in zeta functions. In- 
stead, the splitting-weighted density of states, 

f{E)=J2^EnS{E-En). (1) 



will provide an effective vehicle for evaluating splittings 
directly from complex orbits. Analysis of f{E) in terms 
of periodic orbits follows from approximating it by the 
difference. 



fiE)^N+{E)-N^iE), 



(2) 



between the staircase functions N± (E) for even and odd 
states. The approximation is valid if we use a resolution 
in energy that is much greater than AEn- We next note 
that N±{E) = -(l/7r)Im dE' g±{E'), where g± are 
traces of the symmetry-projected Greens functions M 
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E — En 



(3) 



The traces receive contributions from orbits 7 which ei- 
ther close simply in phase space {g-y = I) or close after 
applying a reflection operation {g-y = Rx)- The group 
characters are x±(/) — 1 and x±{Rx) — ±1, S.y is the 
action and Ay is obtained from the stability of the or- 
bit. Real orbits cannot cross from one well to the other 
and necessarily correspond to g^ = /. Their contribution 
vanishes when we calculate the difference, a fact which 
is true to all orders in h. f{E) is then approximated by 
orbits corresponding to g^ = Rx, which are necessarily 
complex, and from which we can extract the exponen- 
tially small splittings. 

After substituting the contributions of complex orbits 
into the trace formula and integrating, we arrive at the 
following sum , 



fiE) 



Im^ (3y 



iS-y/h 



y/- det(A/^ - I) 



(4) 



Here S-y and are respectively the complex action and 
monodromy matrix of 7. (M^ includes, if necessary, a 
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factor representing linearisation in the surface of section 
of the symmetry operation gj.) We use a complex square 
root in the denominator and the ambiguity in sign is de- 
termined by following the evolution of the square root 
in the complex plane — this mirrors the computation of 
Maslov indices for real orbits . The factor 2 reflects a 
degeneracy in the direction of tunnelling and the dimen- 
sionless factor /3-y, explained below, is 1 for most orbits. 

We explore this approximation for the two-dimensional 
potential 



(5) 



When E < 1, the classical motion is confined to one of 
two symmetric wells, leading to splittings. There is a real 
periodic orbit confined to the x-axis which has a bifur- 
cation at Ec = 0.236 such that it is elliptic for E < Ec 
and inverse-hyperbolic for E > Ec (in which energy range 
phase space is predominantly chaotic) . Tunnelling in this 
system is dominated by a particular set of complex orbits 
that are similarly confined to the cc-axis. Their contribu- 
tions are found by first considering the analogous one- 
dimensional orbits in the potential y(a;,0) = {x^ — 1)'^. 
They are then dressed with the complex monodromy ma- 
trix obtained by embedding them in the other dimension. 
This is useful because a very careful analysis of such or- 
bits is possible; the results of which can later be extended 
to more general orbits. 

The simplest complex orbit starts on the x-axis with 
negative kinetic energy and evolves after an imaginary 
time IT to the symmetry-related point on the other side. 
It has real position and imaginary momentum, giving an 
imaginary action S — iK . The amplitude of this orbit is 
reduced by a factor of 2 on reflection from the inside of 
the barrier. This anomalous reflection cocfhcient derives 
from a calculation in which a Stokes constant is calcu- 
lated on a Stokes line and can be understood in 
simple terms by comparison with an exact calculation 
for the inverted harmonic oscillator [0. As a result, 
the orbit contributes with a factor (3^ = 1/2. The mon- 
odromy matrix can be found by inverting the potential. 
The tunnelling orbit then transforms to a real, unstable 
orbit running along the ridge y = 0, whose monodromy 
matrix M is readily calculated. Reverting to the upright 
potential, the monodromy matrix W is obtained from 
M by multiplying the off-diagonal elements by ±«. This 
leaves eigenvalues unchanged. After including a phase 
factor from reflection inside the barrier, we obtain 



1 e-'^/^ 
ME) = 3 



TT ^- det{W - I) 



(6) 



A final complication derives from the fact that we will 
consider seperately states with even and odd y-parity. 
The orbit lies on the corresponding symmetry axis and 
contributes differently to the two parities as a result. The 



amplitude must then be decomposed according to the 
prescription in . 

For comparison with the theory, we also found the 
quantum eigenvalues corresponding to the potential (|^) 
numerically. We used 100 harmonic oscillator basis states 
in both the x and y directions and diagonalised the re- 
sulting Hamiltonian. We worked to quadruple precision 
which allowed us to calculate splittings as small as about 
10^'^'^. All results shown are for h = 0.01. We also 
used the appropriate bases to isolate the four symmetry 
classes. States even and odd with respect to x generate 
the splittings. The symmetry with respect to y is addi- 
tional and we refer to the corresponding parity classes as 
even and odd. 
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FIG. 1. The dots show the quantum splittings multiplied by 
the Thomas-Fermi density of states and the solid curves show 
the results using fo{E) of Eq. (6). "Even" and "odd" refer to 
the y -parity. 



In Fig. 1 we show as filled dots the exact splittings AE'n 
weighted by the Thomas- Fermi density of states po(-E'n) 
— these are dimensionless numbers expressing the split- 
tings in units of the mean level spacing. The solid curve 
represents the prediction in and agrees well with the 
global trend of the actual splittings. There is a devi- 
ation near E = but this is to be expected because 
the tunnelling orbit approaches marginal stability there. 
Elsewhere the agreement is impressive. In particular, the 
theory works for E < Ec^ yci which range the motion is 
mixed and the real one-dimensional orbit stable. There- 
fore, the theory predicts the mean behaviour whenever 
the tunnelling is dominated by isolated orbits, regardless 
of the character of the classical motion in the wells. 

In addition to the mean behaviour, there is an oscilla- 
tory structure in the splittings. To explain this, we con- 
sider orbits obtained by attaching, to the basic tunnelling 
orbit above, real periodic orbits in the wells on either side 



2 



— always with y = 0. We denote by Sq and Tq the action 
and period of the real primitive orbit segment. Restrict- 
ing the sum to orbits which tunnel only once, we allow 
any number of iterations of the real orbit in the first well 
before tunnelling and in the second well after tunnelling. 
Orbits going from left to right with a total action rSo+iK 
occur with a degeneracy r. On doing the trace integral 
in the wells, we get a contribution to g+{E)— g-{E) with 
an amplitude factor rTo/lih). Summation over these con- 
tributions is sufficient in one-dimensional calculations to 
identify poles in g±{E) and obtain individual splittings 
[^,||. For a consistent approximation of f{E), however, 
it is also necessary to add a term with amplitude it /{iK) 
that arises from integration across the forbidden region 

— a contribution which was not included (or needed) 
in the one-dimensional calculations. The integration in 
energy leading to the staircase functions leaves a denom- 
inator d{rSo + iK)/dE which cancels this combination 
of periods. Including a factor 2 to account for the choice 
of starting well, the result is with = 1. This ar- 
gument can be made precise by considering the problem 
of an infinite square well with a finite square barrier in 
the middle The equivalent calculations are exact 
and free from the ambiguities inherent in comparing ex- 
ponentially small quantities in WKB calculations. 

The monodromy matrix for an orbit with r real han- 
dles is conjugate to Mr — WMq, where Mq is the real 
monodromy matrix of a primitive real orbit segment. Air 
has complex entries. Collectively, these orbits contribute 



2 {rtSo-K)/h 

fosciE) = - Im V , 

^ t[V~det(WMS-I) 



(7) 



to f{E). To determine the branch of the square root, 
we write det(Af — /) = v'A(1 — A~^), where M is a 
surface-of-section matrix with leading eigenvalue A, de- 
fined at each point along the orbit. Initially A = 1 and 
it grows along the real axis if we start with the tun- 
nelling segment. With subsequent evolution in real time, 
A moves into the complex plane and, when the real seg- 
ment is unstable, executes a counterclockwise rotation of 
approximately tt with each period Tq, this becoming more 
exact with every iteration. To this evolution should be 
added a counterclockwise rotation of tt at each turning 
point. Asymptotically, each iteration then leads in the 
unstable case to a phase factor (—i)^ in the amplitude, 
reflecting the Maslov index cr = 3 of the real orbit |pT|| . 

There are analogies between adding fosc{E) to fo{E) 
and adding periodic orbit contributions to the Thomas- 
Fermi density of states. In particular the real phase in 
/osc(-£') introduces oscillatory structure to the splittings. 
To make a comparison with exact results, it is useful to 
divide out the global behaviour apparent in Fig. 1 — the 
result should be a fluctuating function with mean value 
1. We compare in Fig. 2 the Thomas- Fermi-density- 
weighted splittings po{En)AEn/fQ{En) (filled circles) 




FIG. 2. As in the previous figure but with the mean /o(-E) 
scaled out. The dots are the quantum data with dashed lines 
connecting them. The lower solid curve is 1 + fosc{E)/ fo{E) . 
The upper dashed and solid curves are the corresponding re- 
sults averaged with a Gaussian of width w — 0.03. Their scale 
should be read from the right axis. 



with the function 1 -I- fosc{E) / fo{E) (lower solid curve) 
obtained by summing over the first ten repetitions r. 

Let us first discuss the energy regime E < Ec where 
the real orbit is stable. Here summation over r results 
explicitly in delta functions — apparent as the first nar- 
row peaks in Fig. 2. (In order that the quantum data 
not be obscured, side-bands around the peaks have been 
graphically suppressed.) The positions give the mean en- 
ergies of states localised near the orbit in the surround- 
ing island, analogous to a quantisation of stable orbits 
obtained in the usual trace formula ll5|l . The peaks have 



weights of the form (2M/To)e~-^/'* JlOll, where A is a 
function of W and Mq. These weights approximate the 
splittings and are similar to the one-dimensional case 
(for which A — 1). In Table I are listed the positions 
of the peaks and the results of integrating under them. 
These predictions are in reasonable agreement with the 
quantum-mechanical mean levels and splittings respec- 
tively. There are also states which do not correspond to 
peaks in the theoretical curve. We presume that these 
are associated with regions of phase space removed from 
the orbit. It should be noted that in the limit E ^ 0, 
the real orbit suffers an infinite cascade of bifurcations. 
Unfortunately, the lowest even state has an energy near 
the first such bifurcation, so the theory cannot be trusted 
and we do not plot the theoretical curve in this range. 

At energies above E^., the real orbit is unstable and 
phase space dominated by chaos. The theory has smooth 
oscillations and there are no longer individual states as- 
sociated with these orbits. However, the theory does re- 
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TABLE I. Exact and semiclassical energies and splittings. 



produce with some detail a marked periodicity in the 
sphttings. Such osciUations have recently been observed 
experimentally by Wilkinson et al in quantum wells 
and interpreted by them as corresponding to enhanced 
tunnelling in states "scarred" by a real orbit. In our 
formalism, we obtain an explicit quantitative prediction, 
but for averaged tunnelling properties and not (in the 
chaotic regime) for individual states. It is hoped that in- 
clusion of more orbits will yield individual splittings. For 
a quantitative comparison, we compare in Fig. 2 the re- 
sults of averaging quantum-mechanical and semiclassical 
rescaled quantities with a Gaussian of width w — 0.03. 
As discussed before, the theory over-estimates the split- 
tings near E — 1; this can in principle be corrected by 
existing theories for bifurcations jl^. Also, the cascade 
of bifurcations at small energies means that the theory is 
not to be trusted in this range either. Other than these 
effects, it is apparent that we are correct with regard 
to the period, amplitude and phase of the oscillations 
and therefore have a rather detailed understanding of the 
splittings. Notice also that, as far as averaged quantities 
are concerned, there is no qualitative difference between 
the regimes in which the orbit is stable or unstable. 

The problem studied here made use of a class of or- 
bits which are essentially one-dimensional — their two- 
dimensional character being carried exclusively through 
the monodromy matrix. This is an artifact of the reflec- 
tion symmetry about the y-axis and it will be interesting 
to follow these orbits as the symmetry is broken, either 
by adding another term to the potential or by adding 
a magnetic field. This would be the first step in un- 
derstanding how to combine the contributions of many 
complex orbits so as to get individual splittings in the 
chaotic case. Another necessary ingredient is some form 
of symbolic dynamics of the complex orbits so that we 
have a systematic way of searching for orbits; without 
this, it is an extremely difficult task to find long orbits. 
A distinct problem is to determine whether the statistics 
of the splittings conform to any universal distributions. 

The form of analysis developed here may help in un- 
derstanding mesoscopic systems such as Coulomb block- 
ade peaks in which electrons tunnel into a quantum dot, 
and on which to date only statistical analyses have been 
brought to bear For this purpose, we note that 



a very natural extension of the theory developed here 
would allow us to find the width- weighted density of res- 
onances for unbounded problems. The formalism might 
also be used to investigate leakage from or between bil- 
liards with small holes. If the holes are small enough, the 
contribution of leaking orbits will be suppressed diffrac- 
tively so that the widths or splittings will decay alge- 
braically with wavelength, as opposed to the exponential 
decay found in tunnelling. 
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